
log using "C:\Users\Madhav Joshi\Box\2021 Research\Why some countries get old leaders\Research & Politics\R&P log.smcl", replace

hist office_entry_age if _merge==3 & sample==1, normal kdensity

//edit COWcode year cname leader ebirthdate borndate years_inoffice year_in yrborn   eindate startdate enddate office_entry_age lag_v2xel_frefair age_70 ///
//lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65  post_coldwar ///
//africa asia ocenia north_america south_america year_1946_55 year_1956_65 year_1966_75 year_1976_85 year_1986_95 year_1996_05 year_2006_15 if _merge==3 & sample==1
sum office_entry_age lag_v2xel_frefair age_65above lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65  post_coldwar africa ///
asia ocenia north_america lag_v2x_polyarchy south_america if _merge==3 & sample==1
//a model cannot have poliarchy and free & fair election
glm office_entry_age lag_v2xps_party  lag_v2xel_frefair lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode) link(log) ml family(gaussian)
linktest

glm office_entry_age lag_v2xps_party lag_v2x_polyarchy  lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode) link(log) ml family(gaussian)
linktest
glm office_entry_age lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode) link(log) ml family(gaussian)
linktest
glm office_entry_age lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode) link(log) ml family(gaussian)
linktest
//quietly glm office_entry_age lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode) link(log) ml family(gaussian)
//coefplot, drop(_cons) xline(0)
margins, at (lag_v2xps_party=(0(.02)1))
marginsplot



//Replicate the results in linear regression 
//a model cannot have poliarchy and free & fair election
reg office_entry_age lag_v2xps_party  lag_v2xel_frefair lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode) 
collin office_entry_age lag_v2xps_party  lag_v2xel_frefair lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1
linktest 
reg office_entry_age lag_v2xps_party lag_v2x_polyarchy  lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode)  
collin office_entry_age lag_v2xps_party lag_v2x_polyarchy  lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1
linktest
reg office_entry_age lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode)
collin office_entry_age lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1
linktest
reg office_entry_age lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode)  
linktest
collin office_entry_age lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1

//Replicate the result with binary dependent variable age 65 & agbove
logit age_65above lag_v2xps_party lag_v2xel_frefair lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode)  
collin age_65above lag_v2xps_party lag_v2xel_frefair lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1
linktest
logit age_65above  lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1, cluster(COWcode)  
collin age_65above  lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65   if _merge==3 & sample==1
linktest
logit age_65above  lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode)  
collin age_65above  lag_v2xps_party lag_v2x_polyarchy lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1
linktest
logit age_65above  lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode)  
linktest
collin age_65above  lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1
logit age_65above  lag_v2xps_party lag_p_polity2 lag_lnCNTS_TerrGuerWarfare  lag_wdi_gdpgr regular_entry fem_leader lag_wdi_popgr lag_wdi_pop1564 lag_wdi_pop65 post_coldwar africa asia ocenia north_america south_america if _merge==3 & sample==1, cluster(COWcode)  
margins, at (lag_v2xps_party=(0(.02)1))
marginsplot

log close
